Observed
# Load results of multiple pairwise comparisons from 12_Organ_Bias.Rmd
## For all species
if (!params$randomized.mean.expr) {
species <- c("LOC","ELU","DRE")
if (!params$all.nonzero.matrix) {
comparisons.mean.expr.list <-
lapply(setNames(species, species), function(sp) {
readRDS(file.path(params$input.rds.path,
paste(sp, "comparisons.mean.expr.rds", sep="_")))
})
comparisons.log2cv.list <-
lapply(setNames(species, species), function(sp) {
readRDS(file.path(params$input.rds.path,
paste(sp, "comparisons.log2cv.rds", sep="_")))
})
comparisons.resid.var.list <-
lapply(setNames(species, species), function(sp) {
readRDS(file.path(params$input.rds.path,
paste(sp, "comparisons.resid.var.rds", sep="_")))
})
comparisons.expr.var.list <-
lapply(setNames(species, species), function(sp) {
readRDS(file.path(params$input.rds.path,
paste(sp, "comparisons.expr.var.rds", sep="_")))
})
} else {
comparisons.mean.expr.list <-
lapply(setNames(species, species), function(sp) {
readRDS(file.path(params$input.rds.path,
paste(sp, "comparisons.mean.expr.nonzero.rds", sep="_")))
})
comparisons.log2cv.list <-
lapply(setNames(species, species), function(sp) {
readRDS(file.path(params$input.rds.path,
paste(sp, "comparisons.log2cv.nonzero.rds", sep="_")))
})
comparisons.resid.var.list <-
lapply(setNames(species, species), function(sp) {
readRDS(file.path(params$input.rds.path,
paste(sp, "comparisons.resid.var.nonzero.rds", sep="_")))
})
comparisons.expr.var.list <-
lapply(setNames(species, species), function(sp) {
readRDS(file.path(params$input.rds.path,
paste(sp, "comparisons.expr.var.nonzero.rds", sep="_")))
})
}
}
Randomized
# Load results of multiple pairwise comparisons from 12_Organ_Bias.Rmd
## For all species
if (params$randomized.mean.expr) {
species <- c("LOC","ELU","DRE")
if (!params$all.nonzero.matrix) {
comparisons.mean.expr.list <-
lapply(setNames(species, species), function(sp) {
readRDS(file.path(params$input.rds.path,
paste(sp, "comparisons.mean.expr", paste0("rnd", params$seed.run, ".rds"), sep="_")))
})
comparisons.log2cv.list <-
lapply(setNames(species, species), function(sp) {
readRDS(file.path(params$input.rds.path,
paste(sp, "comparisons.log2cv", paste0("rnd", params$seed.run, ".rds"), sep="_")))
})
comparisons.resid.var.list <-
lapply(setNames(species, species), function(sp) {
readRDS(file.path(params$input.rds.path,
paste(sp, "comparisons.resid.var", paste0("rnd", params$seed.run, ".rds"), sep="_")))
})
comparisons.expr.var.list <-
lapply(setNames(species, species), function(sp) {
readRDS(file.path(params$input.rds.path,
paste(sp, "comparisons.expr.var", paste0("rnd", params$seed.run, ".rds"), sep="_")))
})
} else {
comparisons.mean.expr.list <-
lapply(setNames(species, species), function(sp) {
readRDS(file.path(params$input.rds.path,
paste(sp, "comparisons.mean.expr.nonzero", paste0("rnd", params$seed.run, ".rds"), sep="_")))
})
comparisons.log2cv.list <-
lapply(setNames(species, species), function(sp) {
readRDS(file.path(params$input.rds.path,
paste(sp, "comparisons.log2cv.nonzero", paste0("rnd", params$seed.run, ".rds"), sep="_")))
})
comparisons.resid.var.list <-
lapply(setNames(species, species), function(sp) {
readRDS(file.path(params$input.rds.path,
paste(sp, "comparisons.resid.var.nonzero", paste0("rnd", params$seed.run, ".rds"), sep="_")))
})
comparisons.expr.var.list <-
lapply(setNames(species, species), function(sp) {
readRDS(file.path(params$input.rds.path,
paste(sp, "comparisons.expr.var.nonzero", paste0("rnd", params$seed.run, ".rds"), sep="_")))
})
}
}